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Abstract. The Born rule is at the foundation of quantum mechanics and transforms 
our classical way of understanding probabilities by predicting that interference occurs 
between pairs of independent paths of a single object. One consequence of the Born 
rule is that three way (or three paths) quantum interference does not exist. In order to 
test the consistency of the Born rule, we examine detection probabilities in three path 
intereference using an ensemble of spin-1/2 quantum registers in liquid state nuclear 
magnetic resonance (LSNMR) . As a measure of the consistency, we evaluate the ratio 
of three way interference to two way interference. Our experiment bounded the ratio 
to the order of 10 -3 ± 10~ 3 , and hence it is consistent with Born's rule. 
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1. INTRODUCTION 

The Born rule is one of the fundamental postulates of quantum mechanics which states 
that if a quantum mechanical state is described by a wavefunction ijj(r,t), then the 
probability of finding a particle at r in the volume element d 3 r at time t is [1] 

p(r,t)d 3 r = ip*{r,t)^{Y,t)d 3 r = \^(r,t)\ 2 d 3 r. (1) 

The Born rule has well described several experimental results, but no experiments (until 
the recent work by Sinha et al. [21 E]) have been performed to test directly the validity 
of this foundational theory of quantum mechanics. Thus, a deviation from the theory, 
if there is any, would not have been evident. Though quantum mechanics has been a 
very successful theory, it still does not fully satisfy our understanding of the universe. 
Therefore, it is very important to take steps towards experimental verification of the 
Born rule. 

As a direct consequence of the Born rule, interference pattern is produced when 
even a single particle travels through two slits. Quantum interference can be stated as 
a deviation from the classical interpretation of probability for mutually exclusive events 
(e.g. paths, slits, eigenstates and etc.) [21 [3]. For instance, quantum interference of 
two paths 1 and 2 is I 12 = P12 ~ (Pi + P2), where Pj is the probability for a path 
configuration i. Similarly, quantum interference of three paths 1, 2 and 3 can be written 
as I123 = P123 — (P12 + P13 + P23 — Pi — P2 — ?3)- According to Born's rule, three paths 
probability is 

Pl23 = |^l + i>2 + -03 1 2 

= P! + P 2 + P 3 + I12 + Il3 + I 23 

= Pl2 + P23 + Pl3-Pl-P2-P 3 . (2) 

Therefore, the Born rule predicts that I123 = [1]. Here we introduce Po to denote the 
probability of detecting particles when all paths are blocked (ideally zero). Non-zero 
value of Po may rise in the actual implementation due to experimental errors such as 
detector noise. Thus, the measured quantity in the experiment is 

I123 = P123 — P12 — P13 — P23 + Pi + P2 + P3 — Po- (3) 

The three paths experiment tests whether Equation ([3]) is zero by observing the 
probabilities resulting from all possible combinations of independent paths being blocked 
and unblocked, and hence validate the Born rule. Note that this experiment is a 
more precise test for Born's rule than an experiment with two paths in which one 
has to measure the non zero interference pattern and compare it with the theoretical 
prediction [2]. We perform the experiment using an ensemble of spin-l/2s in LSNMR 
quantum computer. This experiment does not only examines the Born rule directly, but 
also demonstrates the capability of LSNMR quantum computing for testing fundamental 
laws of quantum theory. 
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In this letter, we report the results of the three paths experiment that exploits 
NMR QIP for testing the consistency of the Born rule. The remainder of the paper is 
organized as follows. Section 2 explains how to represent the three paths using energy 
eigenstates which can be implemented with NMR. Section 3 describes the experimental 
set up. The results of the experiment and the discussion of possible sources of errors 
follow in Section 4 and Section 5. 

2. ENCODING THREE PATHS IN ENERGY EIGENSTATES 

Here we describe how to translate the triple-slit experiment [21 [3] into a form that 
can be implemented with NMR. Consider a four-level system with energy eigenstates 
|0), |1), |2), |3). One can imagine |1), |2), |3) as wavefunctions of a particle in the 
triple-slit experiment [21 |3] that travels through slit 1, 2 and 3, respectively. We can 
translate the slit experiment into one with three qubits by introducing the notation 

7 = 7l727 3 g {000, 001, 010, 100, 110, 101, 001, 111}, 

where 




if path k is blocked 

1 if path k is unblocked 



and write probability as P 7 for a path configuration labeled by 7 . The |0) state represents 
a slit being closed and is needed for the background probability Po, which is introduced 
earlier for experimental purpose. Now, one can encode one of the eight possible path 
configurations into a state 

l^ 7 > = /3|0> + ]T^#>, ( 4 ) 
k=l v<3 

where f3 is determined from the normalization condition. The table in Figure [1] illustrates 
all possible slit patterns that match with path configuration label 7 , and corresponding 
superposition states l^ 7 ). 



7 


Slit Pattern 


1^) 


000 






001 






010 




^(V2|0} + |2)} 


100 




-L(v^|0) + |l)) 


011 


D 


^(|0) + |2) + |3)) 


101 


IB 


^<[0) + [l) + |3» 


110 


IB 


^(|0) + |1) + |2» 


111 


DD 


^(|1) + |2) + |3)) 



Figure 1: Table of all possible path arrangements that can be formed from three 
independent paths and superposition states that encode each configuration. 
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Moreover, we introduce r to parameterize the evolution between the coherence 
creation and detecting interference, r is directly related to the position of the detector 
with respect to some central position in the triple-slit experiment |2] [3]. Suppose the 
state |^ ul )(V> m | evolves under the Hamiltonian H = E |0)(0| + Ei|l)(l| + E 2 |2)(2| + 
E3 1 3) (3 1 for r. Then the evolved state would have the following form: 

IV> n V)) = E^li>, (5) 

j=l v 

where Aj = Ej — E . Born's rule for probability dictates that 

P 7 (r) = |(^ m (T)>| 2 . (6) 

We can analytically calculate P 7 (t) for all 7 and confirm that Equation ([3]) vanishes for 
all r, E , Ei, E 2 and E 3 if the Born rule holds: 

I(r) = P m (r) - Pno(r) - P 101 (r) - P n(r) + P 100 (r) + P„io(r) + Pooi(r) - P oo(r) 

= - + ^[(cos(A 2 - Ai) + cos(A 3 - Ai) + cos(A 3 - A 2 ))r] - + cos(A 2 - A 1 )r] 

-l[l + cos(A 2 - Ai)r] - 2 -[l + cos(A 2 - Ai)r] + 1 + 1 + 1- = 0. (7) 



3. EXPERIMENTAL SETUP 



An ensemble of a molecule with at least two spin-l/2s in LSNMR suits for the 
experiment. Spin-1/2 nuclei possess superposition of up and down states under static 
magnetic field and act as tiny bar magnets. In LSNMR quantum computing, an 
artibrary single qubit gate is implemented by applying electromagnetic pulses oscillating 
at radio frequency (RF pulses) along the plane perpendicular to the axis of the external 
magnetic field. Two qubit gates are achieved by modulating the evolution under 
" J-coupling" . These single and two qubit gates enable universal control in LSNMR 
quantum computing. A large number of identical molecules precesses around the static 
magnetic field axis and forms a detectable magnetic signal allowing one to make a 
measurement. 

The three paths experiment with LSNMR can be illustrated as a quantum circuit 
in Figure |2j Two qubits are used for encoding three paths, and a third qubit is added 
for read out. As a part of the initial state preparation, we perform "RF selection" [7J E] 
in order to reduce inhomogeneity of RF field strength experienced by our liquid sample. 
Then a pseudo-pure state [5) is prepared by an algorithm proposed in [6]. We use 
magnetic gradient pulses along the z-axis (the direction of static magnetic field) for 
labelling the coherence and decoding it to a pseudo-pure state. The output pseudo-pure 
state is a(X <g> 1 00) (00 1 ) where a is initial spin polarization. Unitary gate U prepares 
\ip lu ), then the state undergoes free evolution for time r. We apply a X n pulse at time ~ 
after U so as to "refocus" [5] unwanted interactions between the third qubit and the two 
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Figure 2: Quantum circuit for the three paths experiment: The third qubit is used to 
read out the information encoded in the first two qubits. The deviation part [5] of initial 
thermal state is denoted as pV Pseudo-pure state (X eg) 1 00) (00 1 ) [5] is prepared by an 
algorithm proposed in [6] with magnetic gradient pulses, and the whole preparation 
procedure is represented by PPS. X n pulse applied to the third qubit at | turns off 
unwanted interactions between the third read-out qubit and the two computation qubits. 
Unitary gate U prepares IV' 111 }; and the state evolves for r into \iP 1u (t)). Applying 
unitary V 1 and measuring in 1 00) basis yield the probability P 7 (t). 



computation qubits. V 1 is an unitary operator that rotates the state along |-0 7 ) to |00). 
Application of V 1 yields the final state whose deviation part [5] is pj = a'^X® 1 00) (00 1 ) , 
where a'^ = aP 1 (r). Measuring the spin magnetization at the end of the experiment 
from the LSNMR spectrum (Figure l4dl) is equivalent to the measurement in 1 00) basis: 

M f = Ti[p)(X® |00)(00|)] 

= «;Tr[(X <g> |00)(00|)(X ® |00)(00|)] = a' r (8) 

Similarly, the initial PPS state magnetization is M* = Tr[a(X <g> |00)(00|)(X g) 
|00)(00|)] = a. Then 

^ = ^ = P(r) (9) 

The three paths experiment was performed in LSNMR on a 700MHz Bruker 
Avance spectrometer at 293X. A three qubit molecule was prepared from a sample 
of selectively labelled 13 C tris(trimethylsilyl)silane-acetylene dissolved in deuterated 
chloroform (Figure [3a|) . Natural Hamiltonian parameters that are relavent for the 
experiment is shown in Figure [30 Two 13 C's are used to carry out the computation 
while spectrum of 1 H is measured. In the experiment, we used Gradient Ascent Pulse 
Engineering (GRAPE) (HI |9] numerical optimization technique to find RF pulse shapes 
that implement unitary evolutions U, V 1 , and the pseudo-pure state preparation with 
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Figure 3: (a) Schematic of a three qubit molecule used for the experiment (not in scale): 
A proton ( 1 H) and two carbons ( 13 C) are used for realizing qubits. (b) A table of natural 
Hamiltonian parameters (Hz), Tl and T2 (s): The diagonal elements give the chemical 
shifts with respect to the transmitter frequencies. The off-diagonal elements are the 
J-coupling constants. 



above 99.95% Hilbert-Schmidt (HS) fidelity defined by 

9 = I I 2 (10) 

We used a 20/is square pulse for refocusing. 

Figure H] illustrates an example of LSNMR spectra of 1 H attained from an 
experiment for measuring P 7 (r). For this particular example, 7 = 111 and r = 0. 



4. RESULTS 

In the experiment, we evaluate the quantity 



|Iho(t)| + |Iioi(t)| + |I ii(t)|' v ; 

where I(r) is the three paths interference, and the denomenator is the sum of magnitudes 
of two paths interferences (e.g. In (r) = Pn (r) - P100O) - Poio(t) + p oooO))- In 
this way, one can assure that the experiment is dealing with quantum phenomenon 
inasmuch as the denomenator should vanish in classical case. Moreover, the calculation 
of such quantity is straight-forward in our experimental set up. As mentioned in the 
previous section, P 7 (r) is obtained by normalizing the magnetization of the final state 
measured after unitary gate V 1 with that of the initial pseudo-pure state. Thus we run 
two experiments consecutively, first to measure the magnetization of the initial pseudo- 
pure state and second to measure the magnetization of the state acquired from the full 
quantum circuit (Figure [2]). These two measurements are separated by 25s (about five 
times larger than Ti) in order for spins to re-thermalize. 

We sample e for various r from 0/is to 1900/is with discritization 5t = 100/is. For 
each r, we repeat the experiment ten times providing 200 data in total. We obtained 
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(c) State after U (d) State after V 111 

Figure 4: X H spectrum is used to gain information about the state of two 13 C's. Vertical 
axis corresponds to magnetization amplitude and horizontal axis is frequency of spin 
precession, (a) is the initial thermal state spectrum of 1 H, and the four peaks are due 
to four possible states of other two qubits. Integrating the rightmost peak at the end of 
the experiment is equivalent to measuring in 1 00) basis. After applying PPS (Figure [2]), 
we obtain a three qubits pseudo-pure state (b). (c) shows the equal superposition state 
|^,iii^^iii I Q f £ WQ carbon Finally, (d) is the spectrum we obtain after applying 



the weighted sample mean (WSM) e = 0.007 ± 0.003. WSM is appropriate for the 
data analysis since the size of standard deviation varies for different r. The random 
error is the standard error of the WSM of e. The results are shown in Figure |5j The 
red dots are the average of ten repetitions of the experiment and the size of the error 
bars indicate standard deviations of the average. The black circles represent simulation 
results. The simulation assumes that the GRAPE and hard pulses designed for the 
experiment are implemented with no error under the effect of T 2 and uses Born's rule 
to extract magnetization signal of the final state. 

5. Analysis of Possible Sources of Error 

In this section, we discuss possible sources of error of the experiment. As shown 
Figure l3bj the difference between the Larmor frequencies of Ci and C2 is an order 
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Figure 5: Experimentally determined e: The experiment is repeated ten times for each 
r, giving 200 data points in total. The red indicates the experimental outcomes and the 
size of the error bars represents standard deviation. We acquired from the experiment 
the weighted sample mean e = 0.007 ±0.003. The blue dots and error bars are obtained 
by repeating the experiment (again, 200 data in total) with the different measurement 
method explained in the following section. The results from the second method are 
e = 0.009 ± 0.003. The black circles are obtained from the simulation which assumes 
the Born rule holds and perfect implementations of the GRAPE and hard pulses under 
the effect of T 2 . The green is the outcome of the simulation that takes systematic errors 
such as distortions in the implementation of GRAPE pulses and RF field fluctuations 
described in the next section into account. The systematic errors well explain the small 
deviation of experimentally determined e from zero; red and blue (experiment) overlap 
with green (systematic errors simulation) for most of the r's, with less than one standard 
deviation away when there is no overlap. 



of magnitude larger than the J-coupling and thus we are well into the weak coupling 
approximation. A simulation of the neglected strong coupling terms shows that on our 
time scale for r, they would contribute to e ~ 10~ 16 and hence negligible. 

Next, there are distortions in the implementation of shaped RF pulses j8]; GRAPE 
pulses seen by the molecule in the LSNMR spectrometer do not exactly match to 
what we desire. There are two components to this deviation; random errors and a 
systematic portion that is primarily caused by limitations with probe circuit design. 
The systematic imperfection can be rectified by placing a pick-up coil at the sample's 
place and closing a feedback loop to iteratively correct the RF pulse shapes [8j [TO] . 
This method improves (yet, still not perfect) the closeness of the actual pulse to the 
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desired pulse. Random fluctuation of RF field is inevitable in the experiment. The RF 
variations for 1 H channel and 13 C channel are found to be 0.7% and 0.2%, respectively. 
The RF selection process mentioned previously is very sensitive to RF field variations 
since it is designed to select a subset of the ensemble of spins at a specific nutation 
frequency. That is to say magnetic signal after the RF selection prior to pseudo-pure 
state preparation between two consecutive experiments can vary [llj. We prepared 
100 pseudo-pure states following the RF selection consecutively and observed 0.95% 
fluctuation of magnetic signal in average with the worst case being about 2%. Hence a 
reference pseudo-pure state taken prior to an actual three paths experiment can be an 
erroneous reference. The fluctuation translates to e ~ 10 -3 ± 10~ 4 . 

The green circles and error bars in Figure [5] show the systematic errors due to 
distortions in the implementation of GRAPE pulses and the random fluctuation of RF 
field. 

We performed another set of experiments with different measurement method. 
Instead of taking two experiments to find probabilities as described earlier, we opened 
the receiver at the end of the pseudo-pure state preparation for a short time, and 
opened the receiver again after V 1 so that the reference magnetization and the 
final magnetization are obtained from a single experiment. We hoped to see some 
improvement in this method by removing slow RF fluctuations between two experiments. 
Nevertheless, there was no significant improvement; we obtained e = 0.009 ±0.003. The 
results from this method is indicated as blue in Figure [5j 

There are other possible sources of error such as transient effect from refocusing 
pulses due to their fast varying amplitude profile and disturbance of static field due 
to gradient pulses used for pseudo-pure state preparation. Furthermore, for the three 
qubit molecule TTMSA in LSNMR, average error per gate is found to be ~ 10~ 3 from 
randomized benchmarking in [TTJ. Such gate error contributes to e ~ 10~ 4 . 

6. CONCLUSION 

The Born rule is one of the fundamental postulates of quantum mechanics and it predicts 
the absence of three way interference. We investigated the three way interference in 
order to test the consistency of the Born rule by performing LSNMR experiment using 
a three qubit molecule TTMSA (Figure l3a|) . The quantum circuit for the experiment 
(Figure [2]) was realized by composing GRAPE and hard pulses. We analyzed the 
quantity e, the ratio of the three paths interference to two paths interference, and 
acquired e ~ 10~ 3 ± 10~ 3 . Some of the major sources of experimental inaccuracy 
are listed in Section |5l and the simulation indicate that the small deviation of the 
experimental outcome from the Born rule is well explained by the systematic errors. 
Potential improvements of the experimental results include ensuring the pulse seen 
by the liquid-state sample is as close as possible to the ideal pulse and reducing RF 
inhomogeneity via enhenced probe design. 

Finally, we conclude that the results of our experiment are consistent with Born's 
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rule. Furthermore, we have demonstrated the capability of LSNMR QIP techniques for 
testing a fundamental postulate of quantum mechanics [121 1131 EH]- 
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